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1 Introduction 



The renormalization group (RG) has proved one of the most important tools in 
refined perturbative analyses. For it has been recognized for a long time that ordi- 
nary loop-wise perturbation expansions of important physical quantities are not only 
restricted to "small" values of the couplings but are often rendered useless by the oc- 
currence of "large" logarithms. RG resummation of these logarithms is then crucial to 
establish a region of validity for perturbative results. 

This is the case in the analysis of vacuum stability (VS) in the Standard Model 
(SM), where the loop-expansion of the effective potential (EP) contains logarithmic 
terms. Only after RG summation of these logarithms may the requirement of vacuum 
stability be turned into bounds on the Higgs mass Again, the discussion of Bjorken 
scaling and its violations in deep inelastic scattering (DIS) is reliable only after RG 
summation of the relevant logarithms yielding in turn high precision tests of QCD and 
one of the most accurate determinations of the strong coupling constant 0. 

To apply the established RG techniques in both cases it is essential that in the 
region of interest (large absolute values of the scalar field in the discussion of VS, large 
momentum transfer for fixed Bjorken variable in DIS) there is only one generic scale 
Ai. Then, using some mass independent renormalization scheme such as MS Ai may 
be tracked by the corresponding renormalization scale /i, as it occurs in the combination 
^log(A1//i^) only. Choosing /i^ = removes the potentially large logarithms from 
the perturbation series. Hence, at this scale the perturbative result is trustworthy for 
"small" values of the couplings and yields the proper boundary condition for the RG 
evolution to finite values of h\og{M. / fi"^) . 

However, there may be many generic scales Adi in the region of interest. For exam- 
ple, in the computation of finite temperature EP or in supersymmetric extensions of 
the SM one encounters this problem [^]. But even in the SM there are largely differing 
effective scales near the tree-level minimum. Although the usual VS analyses of the 
SM were concerned with large absolute values of the scalar field far away from the 
tree minimum it is implicitly assumed that the tree minimum is only slightly shifted 
by quantum corrections. For consistency, one should check this assumption; this is a 
highly non-trivial multi-scale problem. The breakdown of the ordinary RG analysis 
of DIS at small and large is again due to the growing importance of generic scales 
other than the large momentum transfer (for a review see ref. [Q). In both instances 
different potentially large logarithms ^log(A^i//i^) occur in the loop- wise perturbative 
expansion which should be resummed in order to get trustworthy results. But as there 
is only one renormalization scale one cannot trace the various Aii at once and remove 
all the log's from a loop-wise expansion at one particular scale. So, although one still 
has a perfectly good RG equation there is no longer a proper boundary condition to 
RG-evolve from. This problem has been recognized by many authors. 

Sticking to the MS scheme the decoupling theorem was used in ref. [0 to obtain 
some regionwise approximation to leading log's (LL) multi-scale summations. Although 
this is perfectly reasonable, one has to employ "low-energy" parameters, and it is not 
clear how to obtain sensible approximations for these low-energy parameters in terms 
of the basic parameters of the full theory. Alternatively, one of us argued that 
one could still apply the standard MS RG equation to multi-scale problems provided 
"improved" boundary conditions were employed. Although some improved boundary 
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conditions were suggested in some simple cases, no general prescription was given for 
constructing these boundary conditions, and no obvious improved boundary conditions 
were apparent for the subleading log's summation. 

Clearly, one must go beyond the usual mass-independent renormalization schemes 
if multi-scale problems are to be seriously tackled. In the context of the effective 
potential we are aware of two different approaches. In ref. ||^ it was argued that 
one could employ a mass- dependent scheme in which decoupling of heavy modes is 
manifest in the perturbative RG functions. Alternatively, in ref. the usual MS 
scheme was extended to include several renormalization scales Kj. While this seems to 
be an excellent idea, the specific scheme in [0] has two drawbacks. Firstly, the number 
of renormalization points does not necessarily match the number of generic scales in the 
problem at hand, as there is a RG scale associated with each coupling. Secondly, 
when computing multi-scale RG functions to n loops one encounters contributions 
proportional to log"~ (ki/kj) (and lower powers). If some of the log(K j/kj) are "large" 
then even the perturbative RG functions cannot be trusted and used to sum logarithms. 
A similar approach to the one of ref. [T^ was outlined in ref. [jll| though no detailed 
perturbative calculations were performed. 

In this paper we adopt a more systematic approach. Using the freedom of Gnite 
renormalizations we introduce a new "minimal multi-scale subtraction scheme" that 
allows for as many renormalization scales k, i as there are generic scales in the problem. 
Hence, removing all large logarithms at scales = M.i in the new scheme standard 
perturbative boundary conditions become applicable. As in the approach of ref. pO| , 
the multi-loop RG functions in this scheme inevitably depend on the renormalization 
scale ratios KijKy However, within our minimal multi-scale subtraction scheme we 
are able to implement a laxge logarithms resummation on the RG functions them- 
selves. Using these improved RG functions the "partial" RGB's corresponding to the 
renormalization point independence of physical quantities allow us then to resum the 
logarithms for any other choice of scales. 

Much like the SM, the calculation of the effective potential near the tree-level min- 
imum of the broken phase (m^ < 0) in the 0(A^)-symmetric 0^-theory is a two-scale 
problem for 1 < < oo. In our opinion, this is the simplest non-trivial multi-scale 
problem in four dimensions, and so we propose to use this model to demonstrate our 
method. In fact, we are able to analytically perform leading order (LO) and next-to- 
leading order (NLO) multi-scale computations in the 0(A^)-model. Surprisingly, this 
analysis indicates that the assumption that the tree-level is not significantly shifted by 
quantum corrections is only valid for A^ > 4. For 1 < A^ < 4 it appears that there 
might not even be a stable vacuum in the broken phase. 

The outline of the paper is as follows. In section 2 we review the standard MS 
RG approach to LL summations in the single-scale cases A^ = 1 and A^ — > oo. In 
section 3 we motivate the idea of two-scale renormalization and introduce our minimal 
two-scale subtraction scheme. In section 4 we compute the leading order two-scale RG 
functions within our minimal prescription. We use these LO beta functions in section 
5 to compute the LO running parameters, which are then used in section 6 to compute 
the two-scale RG improved potential to leading order. In sections 7 and 8 we determine 
the next-to-leading order contributions to the RG functions and running parameters. 
In section 9 we obtain the NLO effective potential. Section 10 is devoted to a discussion 
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of the special case N = 2. In appendix A we collect the values of various constants 
and in appendix B we discuss some relevant two-loop integrals. 



2 Resumming log's in the effective potential 

Let us consider the massive 0(A^)-symmetric field theory with Lagrangian 
where is an A^-component scalar field. Note the inclusion of the cosmological constant 



A [|T2| which will prove essential in the discussion of the RG later (For a nice discussion 
of this point in the context of curved spactime calculations we refer to JT^). 

We are interested here mainly in the effective potential which arises as the zeroth 
order term in a derivative expansion of the effective action T[ip\ 

V[ip] = J d'x i^-V{^) + ^Z{ip)d^ifd^ip + 0{d^)^ . (2) 

As usual T[ip] is the Legendre transform of the Schwinger functional yV[j]. 

A loop-wise perturbation expansion oi V = J2 (i\2n "\/("-'°°p) W\ |T5[| yields in the 

MS-scheme 

24 2 ^ ' ' 
. (3) 

where 



/i2 2/ ^ '4 1/^' 



Mi = m'' + \\^\ M2 = m^ + IX^^, (4) 

and fi is the renormalization scale. The one-loop contribution to the EP thus contains 
logarithms of the ratios Adi/ to the first power and in general the n-loop contribution 
will be a polynomial of the nth order in these logarithms. (The explicit two-loop result 
has been obtained in [0.) The EP is independent of the renormalization scale /x which 
gives rise to a MS RG equation. 

In view of these logarithms the loop-wise expansion may be trusted only in a region 
in field- and coupling-space where simultaneously 

<1, — -^log^«l, (5) 



(47r)2 ' (47r)2 ^ ^ 

conditions which may hardly be fulfilled e.g. around the tree-level minimum of the 
potential, where in the broken phase A^2 = 0, even with a judicious choice of the scale 
/i. Hence, to obtain a wider range of validity one has to resum the logarithms in the 
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In the two limiting cases A^ = 1 and N ^ oo there is essentially only one relevant 
scale involved, Aii for A^ = 1 and M.2 for A^ oo. Setting the renormalization scale 
/X equal to the relevant scale removes the potentially large logarithms at this scale and 
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we may trust the tree-level EP there. To recover the EP at any other scale we then 
use the MS RG 

We next expand the RG functions in powers of h. As the expansion coefficient Z[(p) 
in (0) does not contain logarithms at the one-loop level no anomalous field-dimension 
arises and it is an easy task to read off the other one-loop coefficients from the result 
. For = 1 we have at one loop 



'^^"(4^' '^"^'""(4^' '^^~2(4^' 1^^-"' 



whereas for ^ oo we find 

_ ^^^^ a _ ^NXm^ _ hNm^ 

which are exact in this limit. 

With the use of the RG functions we next recover the running couplings. Setting 
s = log(/i(s)//i), where jj, is the reference scale, we have for A^ = 1 

\{s) = A(l - 3As)-\ m'{s) = m\l - SAs)-^/^ 

A(.) = A - ^ [(1 - 3As)i/3 - l] (9) 

and for N ^ oo 

A(s) = A(l - lNXs)-\ m\s) = m\l - iNXs)-^ 

A(.)=A + ^[(l-iArA.)-^-l]. (10) 

Imposing the tree-level boundary condition the LL approximation to the effective 
potential at an arbitrary scale fi becomes 

V!^^ (A, m^ ^, A; /.) = + ^m'{s,)^' + A{s,) (11) 

where 

= Of A \2 ^Qg^' *2 = 7777^ log 12 

Higher orders may again be systematically resummed giving rise to the NLL, NNLL, 
... approximations to the effective potential [ll7| . 

As the usual RG may cope with one scale only this approach does not allow a 
systematic resummation in the generic case as we have to deal with two relevant scales, 
at least near the tree-level minimum in the broken phase. Therefore, we have to 
generalize the usual RG approach allowing for as many renormalization scales as there 
are relevant scales in the theory, the task we turn to in the next section. 
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3 Two-scale renormalization 



In the previous section we were able to use the renormahzation scale n arising in MS 
to track one relevant scale and to resum the corresponding logarithms with the MS 
RG. This was sufficient to obtain a trustworthy approximation to the EP for = 1 
and oo. To deal with the general case we shall introduce a new set of parameters 
depending on two renormalization scales /t i, /t 2 which allow us to track the two generic 
scales That is, we consider a finite transformation 

-^Ms = -Fa(A; K 1, K2, /i) 

i^hs = "^^i^„2(A;Ki,K2,/i) 

Ams = A + m^FA(A;Ki,K2,/i) 

V^Ms = V5-FV(A; Ki, K2,/i). (13) 

Here, the MS parameters Ams, ^ms, V^ms at scale /i may be regarded as "bare" 
ones as opposed to the new "renormalized" two-scale subtraction scheme parameters 
A, m^. A, (f. 

Our goal is to construct a transformation ([13D with the following properties: 

i) The effective action F, when expressed in terms of the new parameters, should 
be independent of the MS scale /i. 

ii) When k 1 = k 2 the minimal two-scale subtraction scheme should coincide with 
MS at that scale. 

iii) When A^ = 1 or A^ — > 00 one scale should drop and the two-scale scheme should 
coincide with MS at the remaining scale. 

iv) When = Aii the standard loop-expansion should render a reliable approxi- 
mation to the full EP insofar as ^^^A(fi; 1, K2) is "small". 

In order to find a suitable transformations (|T^) with the desired properties we ffist 
study the associated RG's and RG functions. Having obtained a trusworthy set of RG 
functions we turn them into running couplings and an improved effective potential. 

Our starting point is 

Tms [Ams, m^, Ams, V?Ms; /i] = r[A, m^. A, K 1, K 2] (14) 

from which we derive the two RGE's corresponding to variations of scales Kj, where 
the other scale Kj and the MS parameters are held fixed, in much the same way as the 
MS RG is usually derived. Specializing to the effective potential we obtain 

I,,V = 0. I,. = .,A + .,a,A + ,^„,_i_ + .^^|__,5^^|., (15) 
The two sets of RG functions are defined as usual 



dX dm? dA dip 



for i = 1,2. In general they may be functions not only of A,m^ as are the MS RG 
functions but also of n 2/1^1- 
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Note that property ii) requires that the sum of the two-scale RG functions at 
K 1 = K 2 coincides with the MS RG function at that scale 

^P{ki = K2) +2p.{fil = fi2) = P.,MS (17) 

and property iii) fixes the two sets of RG functions in the single-scale limits. For = 1 
there are no Goldstone bosons. Hence, we have to choose the usual = 1 MS RG 
functions as the first set of RG functions, given to 0{h) by eqn. (|^, and to disregard 
the second set of RG functions so that V2 = K,2d/dn2- For N ^ (X) there are no 
Higgs contributions. Accordingly, in this limit we have to disregard the first set of RG 
functions, so that Vi = n id/dni, and to choose the second set as the large A^ MS RG 
functions, given by eqn. (|). 

Let us come back to the general case. As we want to vary k, 1 and k,2 independently 
we must respect the integrability conditions 

[Kid/dKi,K2d/dK2] = [Vi,V2] = 0. (18) 

An essential feature of a mass-independent renormalization scheme such as MS is that 
the beta functions do not depend on the renormalization scale /i. Unfortunately we 
cannot generalize this to the multi-scale case and demand that the two sets of beta 
functions be independent of k, 2/1^1- The point is that the independence of the RG 
functions from the scales id/dn = 0, is incompatible with the integrabil- 

ity condition eqn. (|18|). However, it is still possible to arrange for one of the two sets 
of RG functions, or in slight generalization for a linear combination of the two sets, to 
be K j-independent. Hence, we assume that 

/3a = if^xPl + XP2, = if^m^Pl + 2(3m^P2 

Pa = iPaPi + 2/5 aP2, = iP^pi + 2l3^P2 (19) 

depend only on A, unlike the RG functions ijS in ([161) . Accordingly, their values in 
a perturbative expansion may be trusted for small A whatever the value of k 2/1^1. Pj 
are real numbers subject to Pi + P2 = 1- The corresponding RG operator 

V = p{D^+P2V2 

d d~d~d~d~d 

= Plf^lT. ^P2l^2T, ^ PXT^ + Pm^^-^ + PATTV - Pv'^T- (20) 

OK I OK, 2 OA Om'^ OA 0(f 

commutes with Kid/dKi. To recover the k 2/'^ i-dependence of Vi we use that eqn. 
(|I^) imphes 

[P,A]=0, (21) 

yielding RG-type equations for the sought-after j/5 . We remark that the final "im- 
proved" potential will have a strong dependence on the p^-parameters. Each pj-choice 
corresponds to a different transformation in eqn. (|T3p which satisfies conditions i), ii) 
and iii). Accordingly, we should decide for which values oi pj the transformation ([T3|) 
"best" meets condition iv). In section 6 we will argue that the appropriate choice is 
Pi = 1 and P2 = 0. That is, the first set of beta functions, which track the Higgs scale, 
are independent of k 2/^ 1. 
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4 LO RG functions 

To determine the j/3 we make a perturbative ansatz 

Note that this is not simply a loop-expansion, since although we expand in h we retain 
all orders in t. Rather, we should view eqn. (^21) as a LL, NLL, ... expansion of the 
two-scale RG functions. Hence, we assume the full k j-dependence of ijS, to enter via t. 
This immediately allows us to rewrite 

The corresponding perturbative decomposition of the RG operators becomes 

with analogous expressions for T),T)^"'\ To determine the respective RG-like equation 
for a given order j/?^"^ we need 



./3i" 


(> 3(«) 


/Si"' 


/ dA 




^~da) 







(25) 



Here, we used the form of eqn. ([TsD implying, in generalization of the single-scale case, 
that iPx, iP^ do not depend on m^, A and j/^A not on A. We can write 

= m^A'^7f^(t), .4"^ = A<^+i5j")(t) (26) 

with analogous but t-independent expressions for the (3^'^\ 
At LO we have a = 6 = and eqn. (pT]) reduces to 

[p(o),pf]=0. (27) 
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The corresponding equations for the various j/?^"-* may be read off from eqn. (^). We 
now solve them in turn. 



iPf^ is determined by 



Inserting the further decomposition p6| ) and taking into account the A-dependence of 
t this equation reduces to 



(p,-pi + ci(°)t)|af = 0. (29) 



Hence, af'^ is independent of t 



.PT = A^af\ (30) 



The equation for ijS^^^l is 



D.P5 + P'?§X'^3-.0f'§-/:l = O (31) 

and reduces to 

{p, -p, + d^'h) + = /?f (32) 

Its solution is best expressed in terms of the function / 

m = p^-p' (33) 



and reads 



^P^S=m'\pi'\ (34) 



where 



6S°) = 5(°)4? and B 



h 



= ^ («^°^A^°H0) - aS^0)/3(°)) . (35) 
The determination of is a bit more involved 

The corresponding reduced ODE then reads 

{p. - Pi + a^'h) ^ + 2/3(M°) = 2/3f ^7^°) (37) 
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and is solved by 



(0) r-2B(0) , 



4 (0) 

m , 



(38) 



where 



-3i 



C(°)aS? and C^'^ 



^(a(^^(0)-af(0)7 



;y(o) 
(0)^n^^,(0) 



25(0) - 



l{0) 



.(0) 
-2j • 



As for iP^^^ the trivial boundary condition (see below) implies 



(39) 



(40) 



In this section we have computed the two-scale LO RG functions for the 0{N)- 
model. The results depend on pj as well as the boundary conditions af'\o), f3l^\o), 
'yi^\o), Sf'\o) which determine the RG functions at t = (ie. ki = ^2)- In fact, at LO 
the boundary conditions are uniquely determined by the single-scale limit conditions 
following from requirements ii) and iii) in section 3 



af\0) = 3, /3r(0) = l, 7r(0) 



?{0) 



(0), 



1 

2' 



5^(0) = 0, 



a 



(0), 



(0) = UN - 1), (3i'\0) = UN - 1), 7f (0) = UN - 1) 



(0)/ 



4^0) 



0. (41) 



3V-' -yj /2 V"; 2' 

The LO RG functions for A and (p are independent of pj, and are given by (some 
relevant constants are given in appendix A) 



i/3a 



(0) 



3\' 



/5f = i(A^-l)A^ = 2/?(°) = 0. 



(42) 



However, the LO RG functions for and A still have a marked dependence on pj. As 
mentioned in the previous section, we are eventually going to adopt the choice pi = 1, 
0. For this choice eqns. 



P2 



and 



(0) 



m A, 



4]) and (|38|) reduce to 



m A 



(43) 



i/3a 



(0) 



\ 2Pr = iN-i) 
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l-3t)-^ + |(l-3t)- 



m 



(44) 

respectively. It is clear that the beta functions possess Landau poles at 3t = 1. Thus, 
these beta functions are only trustworthy for 1 ^ 3t. Returning to the general pj-case, 
the beta functions have a Landau pole at pi—p2 = a^^H. To avoid this pole we require 
Pi — P2 ^ a^^h for pi > p2 and pi — P2 ^ a^^H for pi < p2. The case Pi = P2 = 
appears to be pathological. 
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5 LO running two-scale parameters 

The running parameters in the minimal two-scale subtraction scheme are functions 
of the variables 

h Ki{Si) hX K2 , . 

where Ki are the reference scales. Note that t{si) as given in eqn. (^) is in fact Sj- 
dependent, t(si) = log . The running coupling may be expanded in a series 
in h 

oo + a 

Ks.t) = Y.^^^X^^\s.,t) (46) 

with analogous expansions for m^(sj, t), A(si, t), ^9(5,, t). We now insert these expan- 
sions into eqn. (0) and solve for the LO parameters. 

The equation for the leading order running two-scale coupling is 

^ = A(°)^af . (47) 

As af*^ is constant it is easily integrated 



where the boundary condition is A(sj = 0) = A. 
Turning to the running mass we have to solve 



dm 



2(0) 



dsi 



m 



2(°^AW/3f). (49) 



Pi is given in eqn. (|3^) in terms of the function f{t). As to leading order 



t(s,) = A(°)(s,)(s2-si + -) (50) 



the Sj-dependence of the r.h.s. of eqn. ( ^ ) is quite involved. Its integration yields 



2(0), , 2l^!^Y //^, .... 

m {Si)=m { — - — 1 I — J — 1 , (51) 



with B^^-* = ^-fer'^^foT) cind with the boundary condition rn?{si = 0) = m^. Here, f'^^Hsi 
is the function obtained by inserting eqn. (|50| ) into eqn. (|33D defining f{t) 



^ \ P2-P1 J 

and / = f^^\si = 0). Note that if the two scales coincide we have t = and f^^\si = 
-52) = / = 1- Requirement ii) provides us then with a strong check on the correct 
algebra for the running LO and NLO parameters. 
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We finally determine the running cosmological constant from 



dsi 



I 2(0)n2 (0) 

(m ) 7,^ ^ 



(53) 



With the use of the results (|5TD for rv?^^'' and (^) for 7)"'' we obtain 

,2(0), 



,(0) 



(0) 



+ Lf 



A(o)(s,) ~T 

2{0)/„ \\2 



A(o)(s,) ^ ^ A ^ 



where 



with C(°) 



^1' 



(0) 



^(0) 



-(0) 



^(0) 



25(0) 



25(0) 



(0) 



(54) 



(55) 



(0) 



JO) 



(0) 
'72 



and A(si = 0) = A. 



To LO there is no anomalous field dimension, and so the field parameter (p does 
not run. 

The LO running coupling A(°)(sj) has a Landau pole at A(a^°''si + a'2^S2) = 1 
and clearly our approximation will break down before this pole is reached. If we 
let one of the Sj —>■ —00 (ie. the far IR region) while holding the other fixed the 
coupling will tend to zero as A(°)(sj — > —00) oc (— Sj)^^. Also note that the LO running 
coupling is independent of pj which parameterize the class of finite renormalizations 
under investigation. 

The behaviour of the running mass and cosmological constant is more complicated. 
Consider the combination 




1 + 



[a 



(0) 



+ a2^)KPiS2- P2S1) 
P2-P1 + a^^^t 



(56) 



In the limit investigated ^j^"^ is not generally positive unless pi = or pi = L Of 
course, we thereby assume that t is chosen such as to avoid the beta function poles 
in which case P2 — Pi + a^^^t has the same sign as p2 — Pi- This is disturbing because 
in eqns. (|5lD and (^4]) we are required to take non-integer powers of this quantity. 
Thus, unless = or pi = 1 we are faced with the disquieting possibility of complex 
running and A in a region where the running coupling is very small. Fortunately, we 
will see in the next section that a comparison of our pj-dependent improved potential 
with standard two-loop and next-to-large calculations indicates that = 1 is the 
"natural" choice. 



6 LO RG improved potential 

It is now an easy task to turn the results for the running two-scale parameters into 
a RG improved effective potential. Eqn. (P!3| ) yields the identity 



V"(A, A; K 1, K2) = V{\{si),m'^{si), V9(si), A(si); k i(si), ^2(52)), (57) 
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with Ki{si) defined in (145|) . Next, we assume the vahdity of condition iv) in section 3. 
Hence, if 

Ki{sif = Mi{sj) = m'^isj) + kiX{sj)!f'^{sj), ^1 = ^' ^2 = ^ (58) 

the loop-expansion of the EP should render a reliable approximation to the RHS of 
eqn. (pTf). 

To proceed we have to determine the values of Si fulfilling (^8|). Insertion of the 
Kii^SiY from ( |58D into (^) yields a quite implicit set of equations 

^'"w'"^^- '''' 

Since we are meant to be summing consistently all logarithms we have to solve (|59|) 
iteratively 

oo t a 

^^ = E77^4'^HA,..,.f^) (60) 



a- 



in terms of the LO log's 



4°^ = ^^log^, where ^ = = 0). (61) 

This yields contributions to the sj"^ from both the Sj-dependence of the running two- 
scale parameters and from their own h expansion (^61). For later use we also give the 
NLO term of the result 

.^^(A,..,.^) = l^og^^M^^ where 

Mf\s,) = m2(°^(s,) + A;.A(°)(^.)^'- (62) 
To obtain the corresponding series expansion for the RG improved effective potential 

oo t a 

ViX,..-K,) = J2——V^'^\X,..,,K.) (63) 



we approximate the RHS of eqn. ([57| ) with those terms in the minimal two-scale 
subtraction scheme result for the EP surviving when Kj(sj)^ = A4i{sj). To 0{h) they 
are explicitly given by 



2(47r)2 \^ 4 

We finally insert the running two-scale parameters from (|46| ) into the RHS of (|6^) 
with their arguments Sj coming from (^). Accordingly, an expansion in powers of 
h yields contributions to the V^""^ from both the Sj-dependence of the running two- 
scale parameters and from their own /i-expansion. Keeping only leading order terms 
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we obtain the LO two-scale RG improved effective potential in the minimal two-scale 
subtraction scheme 



24 



(0)^ 



(65) 



Let us next examine its properties. In the single-scale limits = 1 and N ^ oo 
eqn. (^) reduces to eqn. (pJ]) for i = 1 and i = 2, respectively. In the general case 
1 < < oo the rn^ip^- and A-terms in eqn. (^) depend on pj which parameterizes the 
class of finite renormalizations under consideration. Comparison with two-loop and 
next-to-large A^ results will provide us now with a natural value for them. 

We have used a two-scale RG to track the two scales Aii and A4.2. Once the two 
log's have been summed up we can set k i = k 2 = /i, ie. we may write our improved 
potential in standard MS parameters. In this way we can compare the improved 
potential (^) with standard perturbation theory. When now inserting the various 
constants and expanding eqn. (|65D in sf^ up to second order 



l^W(A,m2,^,A;/i) 



+ 



24 



1 + 3A.S°^ + - 



+2{N - 1)A24°)4°) + 



{N 



2 JO)^ 



-X's 



+ 2- 



(0) 



A^ 



1 + Xsl 

{N -l)p2\ .2(0) 



-As 



(0) 



3(^2-^1' 



+ 

4 

m 

T 



A^ 



A^ + 2 



3P2 



A^- 1 



2P2 



P2-P1, 



X2J0) (0) 
A Si S2 



2 JO)^ 



N 



P2-V1, 
-A.f + 



{N -l)p2 \ ^2^(0)2 

3(P2-Pi)/ ^ 



A^- 1 / 
+ — I 1 



2P2 



3 

A^- 1 



+ 



6 



, P2-P1 
N+1- 



x2 (0) (0) 
A Si S2 



2P2 



2 JO)' 



P2-P1, 



X's 



+ A 



(66) 



we see that the 0{s[^'') -terms in eqn. ( |66D agree with the logarithmic terms in the one- 
loop result (^. The quadratic, pj-dependent terms in eqn. (^) should be compared 
with the two-loop MS effective potential [[16| 



-^(2-loop) 



1 - lot 



Ml 



+ (N^ - 1) 



XM2 
24 



, M2 



^.XMiM2 , Ml -^2 , , Ml M2 
+ (N-l) — — — ( 1 -log^ -log^ + log^log- 



iV): 
12 



12 



I{Mi,Mi,Mi 



12- 



12- 



{N-l)^-^IiM2,M2,Mi), (67) 
36 
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where /(x, z) is the general subtracted "sunset" vacuum integral discussed in ap- 
pendix B. 

Note that the sunset integrals do not contribute to the m^-terms. When comparing 
the m^-terms in eqns. (^) and (0) it is easy to see that they only agree for pi = 1 and 
P2 = 0. Comparison of the rn^ip^- and yj^-terms is more tricky due to the non-trivial 
sunset integrals. 

We should decompose these integrals into logarithmic and non-logarithmic parts. 
This is not too difficult for /(A^i, A^i, A^i). Unfortunately, the decomposition of 
/(A^2, -A^2, -^i) is not unique. However, as discussed in appendix B it seems natural 
to adopt the following one 

Hx,y^z) = {y + z - x)log^\og^ + {z + X - y)log^log^ 



+ {x + y 



log — log — 



X y 
+ 2x log + 2y log — 



+ 2z log — + "non-logarithmic" terms. 

/i2 



(68) 

Inserting eqn. (|68D into eqn. (|67D we see that the y^^-term agrees with the one in eqn. 
(|5^). The (/^^-terms agree only if pi = 1 and P2 = 0. 

An alternative check on eqn. (|65|) is provided by the large N limit. By construction 
our improved potential will agree with standard large results. Examining the next- 
to-large A^ result |jl9[ we have found that in the LL approximation the m^-terms in eqn. 
(|65|) and in the next-to-large A^ limit expression only agree if pi = 1 and P2 = 0- To 
compare the m'^ip'^ -terms we have again employed a "natural" decomposition of some 
integrals and once again agreement is achieved for pi = 1 and P2 = 0. We remark that 
no other choice of pj may be obtained by simply adopting a different decomposition 
of the relevant integrals. We have been unable to check the 93^-terms since we do 
not know whether it is possible to perform a "natural" decomposition of some of the 
contributing integrals. 

Thus, a comparison of our improved potential with the standard two-loop and 
next-to-large A^ potentials strongly indicates that pi = 1 and p2 = is the appropriate 
choice. This is gratifying, since for this choice one does not encounter the complex 
running parameters mentioned in the previous section. Let us finally write down ex- 
plicitly the two-scale improved potential in the two-scale minimal subtraction scheme 
for this choice of pj 

o^J0) AT-l 



+ 



24 



1 - 3Asi' 



1 - 3Xs'; 



(0) 



3 

A^ 



-As 



(0) 



1 -3t 



2 N-1 
' 3 N+S 



m 
2A 



1 - 3Asi' 



(0) 



A^ 



-As 



(0) 



+ 2 



N -A~ 



3t) 



Af+8 



As 



(0)^ 



As 



(0)' 



1 -3t 



' JV+8 



3t 



- 1 



4 Af-l 
' 3 iV+8 



+ A. 



(69) 
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For t = this result has already been obtained in a different way in ref. §]. In the 
broken phase (m^ < 0) the tree-level minimum is at A^2 = or s^2^ — > —oc. Hence, as 
we approach it log(A^2/-^i) will become large. If we are prepared to trust eqn. (|69|) 
even in the extreme case of the tree minimum itself an intriguing property emerges. 

As long as > 4 the ip^-and m^</9^-terms vanish and the m^-term converges to 
a finite value. As the slope ^^wi^'^^ ~^ —oo) \ the EP takes its minimum in the 

broken phase at the tree-level value and becomes complex for even smaller y)^- values. 
But for 1 < A^ < 4 the m'^-term, and hence V"*^*'-', diverges to minus infinity. This 
indicates that for these values of A^ there is no stable vacuum in the broken phase. 
Note especially that for A^ = 4, ie. the SM scalar boson content, the divergence is 
softer but still there, as the penultimate term in eqn. (pQf) becomes a logarithm 



V^'' = "1^(1 -3^)^ log (l-^l+A. (70) 



7 NLO RG functions 

The LO results of the last two sections have already been obtained in a less general 
form in ref. P] based on the use of the MS RG and the conjecture that the correct 
boundary condition at /i^ = Ai2 are given by the A^ = 1 result (|TTp. But using those 
techniques it appeared to be impossible to go beyond LO. The finite renormalization 
(|l^), introducing the appropriate number of renormalization scales and the correspond- 
ing RG equations (0), allows us to overcome this problem in a systematic manner. To 
show the strength of this technique we now determine the NLO RG functions and in 
the next section the corresponding NLO running parameters. 
To NLO eqn. (|T]) yields 

] + =0. (71) 

The corresponding equations for the various are obtained with the use of eqn. 
(p5|). We now solve them in turn, 
is determined by 

CM" + /3f - + - ,Pf'§-,P? = 0. (72) 



Proceeding in an analogous way as in obtaining the LO RG functions in section 4 we 
easily obtain the solution 



.tf ^ = (73) 



where 

«r^w=aS?+4-v^^(t); 

4P=i«4? and i« 



a 



4-^ = ^(«^V^(0)-af(0)««). (74) 
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The equation for is 

n + Ml) _ .f)W_^m 3(1)_^.^(0) _ .^(0)_^3(1) _ n (^r:) 

with the solution 

= ^ (76) 



where 



The equation for i/3^l^ becomes quite involved 



(1) 3(0) d (1) (1) 9 ~(o) 3(0) 5 ^(1) ^(1) d ~(o) 



+Pa 9;^»^a iPx q^Pa ^Prn^Q^2'P^ ''^"'^ dm^^ " ^ ^ 



After some algebra we find the result 



/3?)=m^A7r, (79) 



where 



c 



ii) = (7«al°) and C« = ^ 



(1) _ 2^(0) (1) (1) _ _ 7(1) (0) 



c^:^ = 7i^^(0) - eg) - c^^) - ci^) - c« = -2A(^)B(o)c£). (80) 

To NLO the anomalous dimension is non-trivial and we have to determine if3^^^ 
from 

The solution is easily obtained 

iPS^-X'Si'\ (82) 
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where 



4) = D«aS? and D« ^^'^ 



4^^ = ^(«^°^5f^(0)-af(0)5«). (83) 



So far we have not specified the values of the NLO boundary constants a-^''(0), 
/5f ^(0), tI^^O) and 5f ^(0). At LO the relevant constants were completely determined by 
the single-scale limit conditions following from requirements ii) and iii). Unfortunately 
they do not anymore uniquely fix the NLO constants. For suppose we expand a[^''(0), 
A^^^(O)? cLnd Sl^\o) in powers of (A^ — 1). Then the large limit condition 

forbids any terms proportional to (A^ — 1)^ and higher powers of (A^ — 1) 0], and the 
A^ = 1 limit condition fixes the contributions proportional to (A^ — 1)°. However, these 
limits tell us nothing about NLO terms proportional to (A^ ~ !)• Of course, we still 
have the condition that the sums of the two sets of RG functions at t = are just the 
usual MS RG functions, ie. if3^^\t = 0) + 2p^^\t = 0) = [3%°"^ In MS /^^'J''^ = 
and the other two- loop beta functions can be found eg. in ref. |TB[. Putting all this 
together we have 

a«(0) = -ll-[l + q,]^N-l), a«(0) = gi(A^-l), 
/3i'\0) = -i-[^ + g2](Ar-l), Pi'\0) = q^iN - 11 



7^(0) = g3(A^-i), 7r(0) = -g3(A^-i), 

-+\- 

12 ~ 1-36 



= T^+[4 + g4](A^-l), Si'\0) = -q,{N-l), (84) 



where qj are real numbers which are independent of A^. We shall comment further on 
sensible choices for qj in the discussion of the NLO effective potential in section 9. 



8 NLO running two-scale parameters 

Using the LO results and the set of RG functions obtained in the last section we 
now calculate the NLO running two-scale parameters, which will be used to construct 
the NLO effective potential. 

The equation for the next-to-leading order running two-scale coupling is 

^ = 2AW«f A« + A(°)'ap\ (85) 
dsi ' 

With the use of the results (|i8| ) for A^°^ and ( [73D for a^^^ we may integrate this equation 
and find 

(1) , (1) 

Abnvp 4(1) - "i +"2 
Above A - ^(0)^^,0) • 
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Turning to the NLO running mass we have to solve 



dm 



2(1) 



dsi 



dX 



The integration of this equation is quite involved and yields 



+ M 



(1) 



A(°)(s,) A 



/(0)(, 



M^'hog 



+ M'i^ log 



'AW(£0 



f 



X 



where 



Ml 



(1) 



M. 



(1) 

2 

r(i) 



(1) _ 



5(1) _ 5(0)^(1)^ 

5(1) _ 5(0)^(1) _ ) + A^^\m - log/, 



Above = 



(o)_L^(o) ■ 



The NLO running cosmological constant is determined by 



(0) 



dsi "' " ' ^ " ^ ^ 9A 

With the use of the various results above we obtain after a tedious computation 



A«(sO = XL[ 



(1) 



AL?) 



A(o)(s,) ~T 
(!!!!!!M!f(o)^, 



' X ■' 



+ 4^^ 



/(0)(, 



m 



f 



25(0) 



L^e^ log 



7(0) 



/ 



, ^(i)K!!m)!w 



A 

'x^ 
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where 



(1) 



r(l) 
-^3 



'(1) 



1 

T 



2Mj 



(1) 



(7(1) 



r(l) 
-^2 



/ 



-(0) 

^2 ) 



-(1) 



2^(0) 



-M, 



(1) 



25(0) 

C'(o) 



C(i) 
25(0) 



25(0) 

+ i(i) f ^(0) - C(o) 



- f 25(1) _ 25(0)^(1) + 
1 V 25(0) 



5(0) 



2^(0) 



'(1) 



(1) 



25(0) _ 1 
^(0)^(1)- 41), 41) = C'(°)A(i) 



(5(°)-5(0)) log/, 



2^(0) 
25(0) _ i 



-M, 



(1) 
5 ' 



(1) 



-C(°)A(i) + l(^) - L 



(1) 



(92) 



Above C(i) 



— foj- We remark that most individual integrals occurring in the 

'-'1 +°^2 

computation of not only A(i) but also A(o) and m^*^^"* yield hypergeometric functions 
and that only the respective sums of those are again expressible in terms of elementary 
functions as given above. 

Finally we determine the non-trivial NLO running of ip{si) 



dsi 



V°)A(°)\(^ ^{s. = 0) = ^. 



(93) 



The integration of this equation is straightforward and yields 

(^(i)(.,) = -^dW\x(o)^s,)-\ 



+ ¥^(5(1)- 5(1)) 



A(°)(g, 
/(°)(3. 



(94) 



where D^^^ 



4°) +4°) 



It is easy to see that X^^\ m^(^) and ^9(1) vanish for > 1 in the limit of one 
Si ^ —00 while holding the other fixed. A(i) will tend to a finite value in this limit 
only for A^ > 4. However, it will diverge for l<A^<4ifpi = l and S2 —00 with 
the same rate as A(o) due to the first two terms in (pi]). 



9 NLO RG improved potential 

It is straightforward to extract the two-scale NLO potential from the standard per- 
turbative boundary condition eqn. ( |64D 



24 



6 
20 



1 2(1)/ (0)n 2 I 2(0)/ (0)x (1)/ (0)x , .(1)/ (0)x 



+ 

+ i:4°'(^r)S + ./3S(^f)^ + .<'(^f 



2 r ,^4 ,„2 

(l)/\ . J0)^ 



24 ' ™^ ^ ' 2 



^r(A,...;.r) 



+ ' Y ' +{N-i y ' Y ' |. (95) 

The different contributions come from the expansion of the running two-scale parame- 
ters, from the expansion of their Sj-dependence and from the exphcit one-loop term in 
(0). In practice, we immediately set pi = 1 and p2 = as has been done in the LO 
result. 

Next, we fix the values of Qj used to parameterize the NLO boundary functions in 
eqn. ( p^ by comparing the g^-dependent NLO potential and the NLO Z (ip)^^^ -function 
with the corresponding standard two-loop results. This immediately fixes qs = and 
hence 7i^'*(0) = 0. The value of ^4 depends on how we decompose the two-loop integral 
J for Z(</))*^^"'°°p) into its logarithmic and non-logarithmic pieces. 

In order to determine the "natural" decomposition of this integral it is helpful to 
consider the general integral J{x, y, z) as given in appendix B. It is symmetric in x, ?/, z. 
Accordingly, a natural decomposition should respect this property. In fact, there is only 
one decomposition which does this 

X y z 

J(x, y, z) cx log — + log — + log — + "non-logarithmic" terms. (96) 
/i^ ji^ 

We are interested in the case J(A^2, -A^2, A^i) and so we choose the coefficient of the 
log (A^ 2 //i^) -term in J {M2, M2, Mi) to be twice that of the log(A^i//i^)-term. This 
implies that the coefficient of (A^ — 1) in 52^'' (0) must be twice the coefficient of (A^ — 1) 
in or g4 = 

To determine qi and q2 we need the subleading logarithms in I{jvi2, M.i)- 
Using the decomposition ( p7D yields qi = and q2 = — ^. Putting this all together, 
the complete set of boundary functions are 

«!'^(0) = -f-i(Ar-l), 4')(0) = -i(Ar-l), 

Pi'\0) = -i-^(Ar-l), /?W(0) = -^(A^-1), 
XI) /n^ _ n .,(1)^ 



7^(0) = 0, 7^(0) = 

T2~^ik^^~^^' — 54 



^^^(0) = 3^ + T^(iV-l), 4'\0) = MN-l). (97) 



The behaviour of the NLO contribution is of most interest around the broken phase 
tree-level minimum, where M.2 = or 82'^ — > —00. As in the LO case all the terms in 
eqn. (pSD will vanish or converge to a finite limit if A^ > 4. But for 1 < A^ < 4 A^^) and 
2/9^"^ ■ ^2^^ will diverge. It is easy to check that they diverge at the same rate as 

A(o) 

m 

the LO analysis. However, as the NLO divergence is suppressed by a factor <ti 1 
qualitatively nothing will change. 
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10 The relevance of = 2 



From diagrammatic considerations we would expect the m^-term in the RG improved 
potential to have a certain exchange symmetry in the N = 2 case. For the case N = 2 
the m^-contributions are invariant under the exchange of Higgs- and Goldstone-lines. 
We would therefore expect that for k i = k 2 = yU the m^-terms in eqn. (|69D should be 
symmetric in s^i^ and sf\ A glance at eqn. (|69|) in this case, 



m 



(1 -3Asi' 



(0) 



2A 

+ other terms, 



;i - f )- 



^ +2(1 



10 



As 



(0)^ 



-3 



(9^ 



clearly shows that the m^-term is not symmetric in s^^^ and s'^\ We find it somewhat 
disturbing that our approximation scheme does not respect this symmetry. 

We know from section 6 that eqn. (pHD matches standard perturbation theory 
through to two loops. Therefore, this sf'^ ^ 82'^ symmetry must go down beyond the 



two-loop level. Expanding eqn. ( ]98[ ) in powers of Si'^ and S2'' up to 0{X^) 



JO) 
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2 JO) fo) 



joy 
■^2 



+A^ 



{OYM . (o)(o)^ . 5 JO) 



+ s 
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:^3f 10j0)4 , 20J0)3 (0) 4 JO) 2 (0)2 20 JO) (0)3 10 JO)' 

,^4{ 22 J0)5 50 JO) ^ (0) 80 JO)^ (0)2 
-|-/\ I 3 "T g i'2 -|- 2Y^1 ^2 



20 joy 



+ 9 ■^l 



178,(0)^(0)^ 



+ 27 ■^l 



I 986 „(0)^ 
+ 135 "^2 



we see that the Si'^ ^ S2'' symmetry survives at three and four loops, but breaks 
down at Gve loops. So we see that the failure of our approximation to observe it only 
appears at quite a high order in perturbation theory. We are unable to explain this 
phenomenon further. 



(99) 



11 Conclusions 

In order to deal systematically with the two-scale problem arising in the analysis 
of the effective potential in the 0(iV)-symmetric </)^-theory we have introduced a gen- 
eralization of MS. At each order in a MS loop-expansion we have performed a finite 
renormalization to switch over to a new "minimal two-scale subtraction scheme" al- 
lowing for two renormalization scales k, j corresponding to the two generic scales in the 
problem. The MS RG functions and MS RGE then split into two minimal two-scale 
subtraction scheme "partial" RG functions and two "partial" RGB's. The respective 
integrability condition inevitably imposes a dependence of the partial RG functions on 
the renormalization scale ratio k, 2/1^1- Supplementing the integrability with an ap- 
propriate subsidiary condition we have been able to determine this dependence to all 
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orders in the scale ratio and have obtained a trustworthy set of LO and NLO two-scale 
subtraction scheme RG functions. With the use of the two "partial" RGE's we have 
then turned those into LO and NLO running two-scale parameters exhibiting features 
similar to the MS couplings such as a Landau pole now in both scaling channels. Using 
standard perturbative boundary conditions, which become applicable in the minimal 
two-scale subtraction scheme, we have calculated the effective potential in this scheme 
to LO and NLO. To fix the remaining renormalization freedom we have compared our 
results with two-loop and next-to-large N limit MS calculations. As a main result we 
have found in both LO and NLO that for 1 < A?" < 4 there is no stable vacuum in the 
broken phase. 

The vacuum instability in the broken phase of the 0(A^)-modcl raises immediately 
the possibility of a similar outcome in a multi-scale analysis of the SM effective po- 
tential. As the method outlined generalizes naturally to problems with more than two 
scales we are in a position to investigate systematically the different possible scenarios. 
Before turning to the SM itself it proves useful thereby to study the effects of adding 
either fermions as in a Yukawa-type model or gauging the simplest case of A^ = 2 as in 
the Abelian-Higgs model. The Yukawa case will either be a two- or three-scale problem, 
depending on whether one includes Goldstone bosons or not. The Abelian-Higgs model 
in the Landau gauge will be a three-scale problem to which the methods in this paper 
are easily extended. Now one has three integrability conditions [Di, Vj] = and one 
must impose three independent subsidiary conditions analogous to [/t i, Pi] = 

which we used in our 0(A^)-model analysis. Note that for the general n-scale problem 
one would have — 1) integrability conditions which should be supplemented by 
— 1) subsidiary conditions. The question of whether fermions or gauge fields 
may stabilize the effective potential for small A^ in a full multi-scale analysis is under 
investigation. 

We do not see any fundamental problem in applying the framework presented here 
to multi-scale computations at finite temperature, to the analysis of the multi-scale 
EP in supersymmetric extensions of the SM or to a full multi-scale treatment of DIS 
problems in the regions of very large or small xb- The necessary adaptions are under 
investigation. 
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A Values of various constants 

Here, we give the values of various constants appearing in the paper. We quote them 
for the choice pi = 1 and p2 = 0. 



N + 2 

N + 8' 



2{N + 8) 



3N 



(100) 
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I, C(o) = ^ (101) 

_3iV+14 ^ 5{N + 2) 

N + 8 ' 6(iV + 8)' 

0, ^ + \ (102) 

12(Ar + 8) ^ ^ 

17(iV + 8) ^(,) 5(iV + 8) 



81 ' 162 
V-M 
324 



0, ^(1)^^^ (103) 



19(A^ + 8) 
486 ' 

{N - 1) (19iV2 - 578N - 2600) ^ {N - 1){3AN'^ + + 2178) ^ 



486(iV + 8)2 243(A^ + 8)^ 

(iV-l)(17iV + 46) ^(1)^ 3iV + 14 



243(Ar + 8) ' * 3(Ar + 8)' 

2(Ar- l)(3Ar + 14) 
3(iV + 8)2 

19(Ar + 8) _ (TV - l)(19iV=^ - 578iV - 2600) 

486 486(Ar + 8)2/ 

{N - l)(34iV2 + 544iV + 2178) 



(105) 



243(A^ + 8)2/ 



log/, 



38(A^ - l)(iV + 8) ^ 2{N - 1)2(19A^2 _ 573^ _ 26OO) 



243(A^ - 4) 243(A^ - 4) (A^ + 8)^/ 

4(A^ - l)2(34Ar2 + 5447V + 2178) (1) _ 35(7V + 8) 

243(7V-4)(Ar + 8)2/ 3 - ^ , 

{N - 1)(19A^2 _ ^YSN - 2600) {N - 1){3AN^ + 544Ar + 2178) 
486(iV + 8)2 243(iV + 8)2 

(N -1)(N^- A2N^ - 360N - 760) 34(7V - 1) , , 

H -TT, -^Ogf, 



log/, 



9{N + 2){N + 8y 81 
(iV- l)(17A^ + 46) + 

486(Ar + 8) ' ^ ~ ~6(iV + 8)' 
2(jV-l)(37V + l4) _ 2(iV-l)(3Ar+14) 

3(7V + 8)2 ' - (TV + 8)2 
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B The integrals / and J 

Here, we list some useful formulae regarding the two-loop integrals / and J. The 
general unsubtracted scalar sunset integral in D dimensions is defined as 



lDix,y,z) = J 



Di 



(107) 



A full calculation of this integral is rather involved [^. However, there is a formula 
in ref. |^T|] which nicely splits the integral into a very simple, for D = A divergent 
expression plus a finite term which is proportional to Id~2{x, y, z), ie. the same integral 
in two lower dimensions. 



lD{x,y,z) 



(47r 



{D-2)iD-3) 



{x-y- z){yz)2 



+ {y — z — x){zx)'2^ + {z — X — y){xy)' 



D-2 



- {A7r)-\x^ + y^ + z^ -2xy -2yz -2zx) lD-2{x,y,z). (108) 

Since the last term is finite we regard it as a "non-logarithmic" term and ascribe the 
logarithmic terms purely to the simple, divergent piece. The renormalized I{x,y,z) 
referred to in the text is then given as 



I{x,y,z) 



FP 



(47re-V)'^(/4-2.(x,y,, 
1 



(47r)2e 



{K4-2e{x) + K4-2e{y) + Ki_2e{z 



where FP denotes the finite part, 7 is Euler's constant and 

d'^k 1 



Knx) 



(27r)^P + 



(109) 



(110) 



The KD-teims in eqn. (B.3) are due to the subtraction of one-loop sub-divergences. 
The unsubtracted Jr){x,y,z) is defined as 



JD{x,y,z) 



d 



d^k d^^l 



Di 



1 



dp^ J (27r)^ (27r)« (P + x){P + y){{k + I + pf + z) 
The renormalized J{x,y,z) which enters into Z {lp)^^'^°°'^^ is simply 



p2=0 



J{x,y,z) 



FP 



(47re-V)'^^4-2.(x,i/,z) 



:ill) 



(112) 



Above, the x, y, z are the (masses)^ on the three internal lines. 
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